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We apply the dual-varibles approach to the problem of the optical response of an disordered film of 
metal particles with dipole-dipole interaction. Long range dipole-dipole interaction makes the effect 
of spatial correlations significant, so that dual- variables technique provides a desirable improvement 
of the coherent-potential results. It is shown that the effect of nonlocality is more pronounced 
for a medium-range concentration of the particles. The result is compared with the non-local 
cluster approach. It is shown that short-range correlations accounted in the cluster method reveal 
themselves in the spectral properties of the response, whereas long-range phenomena kept in the 
dual technique are more pronounced in the k-dependence of the Green's function. 



I. INTRODUCTION 

Disordered systems of a various physical nature with 
different typical lengthscales play an important role in 
condensed matter physics. Among the systems with a 
disordered atomic configuration we find various alloys, 
amorphous semiconductors and diluted magnetics. Dis- 
ordering at a nanoscale level is typical for arrays of semi- 
conductor quantum dots and metallic heterostructures. 
Depending on the lengthscale, the disorder reveals it- 
self in different experimental observables. For systems 
with an atomic length scale disorder, the electronic prop- 
erties are the most interesting. For random arrays of 
nanometer-size and larger particles, the study of their 
optical response is a remarkably active field of research. 

There is a number of common approaches to describe 
the different phenomena in disordered systems. In par- 
ticular, the approach of the coherent potential approxi- 
mation (CPA) has been used for the description of the 
electronic system in the substitutional alloy modeP^, but 
it has also been applied in the description of the opti- 
cal properties of disordered nanostructures. Neverthe- 
less, there is a considerable difference in the type of the 
interaction (short-range hopping for electrons and long- 
range Coulomb interaction in the optics of nanostruc- 
tures). One should also note that optics offers much 
larger set of observables; for example, the generation of 
high optical harmonics allows experimentalists to study 
many-particle Green's functions of photons. For a the- 
oretician working in the field of electronic structure, a 
description of the optical response provides an additional 
possibility to validate existing calculation schemes and it 
can be a good playground for the development of new 
methods. 

In this paper we consider a model describing disor- 
dered films of metallic particles. The optical properties 
of small metal particles in a dielectric host have already 
been investigated for a very long time. Inhomogeneity 
is common for these materials - the shape, size and po- 
sitions of inclusions in the matrix can be distributed in 
a random way. There exists a sizable number of mod- 
els and theories devoted to the calculation of the optical 
response of such metal-dielectric composites. Nowadays, 



the problem of theoretically describing and calculating 
the optical properties of such systems have gained inter- 
est because of progress in the technology of preparing 
novel materials^. The possibilities of modern technology 
to fabricate samples with well controlled parameters al- 
low to perform an accurate experimental verification of 
theories in this field. 

The solution to the problem of light scattering by a 
single spherical particle in a host material was provided 
by Mie in the beginning of the 20th centurjffSl^ Metal par- 
ticles exhibit a well-known Mie resonance, corresponding 
to excitation of the collective mode of the conduction 
electrons. For ensembles of metal spheres, the Mie reso- 
nance peak is broadened and shifted due to the electro- 
magnetic interaction of particles and disorder effects. In 
the three dimensional case, these effects lead to certain 
spectral behavior corresponding to an effective dielec- 
tric constant of the inhomogeneus medium. This result 
can be achieved by the so-called mixing approximations, 
the most well-known were developed by Bruggeman'^ and 
by Maxwell GarnetlP. In Bruggeman's method, all con- 
stituents of the composite, including the host material, 
are considered as spherical inclusions in a homogeneous 
effective medium. The Maxwell Garnet theory (MGT) 
considers the effective medium for the inclusions, but not 
for the matrix material. In both approaches the dielec- 
tric constant of the effective medium is obtained self- 
consistently. Through a number of more sophisticated 
effective-medium approaches, the CPA can be shown to 
be the best local theory, as it accounts of all local dia- 
grams. For the planar systems such as metal island films, 
the effective dielectric constant makes little sense. How- 
ever, the concept of an effective medium is still fruitful, 
and methods of the CPA family are widely accepted. 

A particular case of a system where the CPA is easy 
to formulate is the case of a periodic lattice with metallic 
particles randomly placed at node^ (the so-called lat- 
tice gas). The CPA corresponds to a substitution of the 
partially filled lattice with a regular one, having particles 
with a self-consistent permittivity at each nod^^SllillH 
Physically this means that fluctuations of electromag- 
netic field at different sites arc assumed as uncorrelated 
(that is, fluctuations are local). Usually, the Drude model 
is used for the frequency dependence of the permittivity 
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of the particles, and a dipole-dipole interaction between 
the particles is assumed. 

There are a few arguments for physics beyond the 
single-site CPA. Firstly, in some systems there are corre- 
lations due to grouping of the disordered particles. For 
example, a fractal nature of clusters is assumed in a series 
of workJSlllljllllill, The idea that in a set of randomly dis- 
tributed particles some specific configurations of particles 
placed close to each other gives a significant contribution 
to the total polarizability is applied irfSl. Another case 
are systems in which the particles are distributed inde- 
pendently, but spatial correlations appear in the electro- 
magnetic field at different particles. This is particularly 
the case for systems with a high filling factor, whereas 
the concept of local correlations is more suitable for low 
filling factors. 

For the description of non-local field correlations, there 
are strong insights from the theory of electronic struc- 
ture in the Anderson model with (5-correlated disorder, 
since the mathematical formulations of both problems 
almost coincide. While the CPA was successfully used in 
a number of calculations for the electronic properties of 
metallic alloys ab initio, certain peculiarities in the elec- 
tronic density of states are attributed to non-local cor- 
relations. A trivial generalization of the single-site CPA 
is to introduce a "supercell" containing several atoms, 
and to suppose that correlations are located in such a 
real-space cluster. However, such an approach clearly 
violates the translational invariance of the lattice and it 
also does not provide a significant improvement to single- 
site theorjI^llEl. It is much better to formulate a theory 
operating in quasimomentum space, as is done in the dy- 
namical cluster approximation (DC A). Roughly, in the 
DCA the Brillouin zone is divided to several parts, and 
the self-energy is supposed to be k-independent within 
each of these parts. Self-consistent calculations for cer- 
tain clusters are used to estimate the self-energy; the 
number of atoms in the cluster equals the number of 
parts of the Brillouin zone. If this number is equal to 
one, the method passes into the CPA. Originally the DCA 
was developed for ordered correlated electronic systems 
(such as the Hubbard model) to include non-local correc- 
tions to the dynamical mean-field approximatiorf^. The 
idea of the DCA was modified and applied to the An- 
derson model of disordered systema^^lj and it is known 
as the non-local coherent potential approximation (NL- 
CPA). It is shown that this method can be reformulated 
as a unique reciprocal-space theory of disordei^^. 

Numerical simulations of the optical properties of dis- 
ordered metal films were also performed, and the results 
of these simulations are useful to test our theoretical in- 
vestigation, because of the possibility to test the systems 
with defined parameters. Using Monte-Carlo computer 
simulations, the dielectric properties of hard spheres of 
point-like dipoles fluids are investigated and deviation 
from the Clausius-Mossotti-Lorenz-Lorentz (CMLL) for- 
mula is investigated irP, and the value of the correction 
parameteiffS has been calculated. For the same type of 



system with different concentrations the effective sus- 
ceptibility is calculated in the worl^I^, the results are 
compared with the results of the MGT approach. IiP 
a Monte-Carlo method is applied for the calculation of 
the response of disordered lattice systems. A spectral 
approach for the calculation of the susceptibility of the 
two-components composites is used in^. 

Apart from the remarkable progress in theoretical 
methods as well as in experimental techniques, a rela- 
tively small number of theoretical predictions have been 
quantitatively compared with experimental data so far. 
This is because of the difficulty in preparing a well- 
controlled sample exactly corresponding to the system 
under theoretical consideration (a lattice gas, for in- 
stance). On the other hand, the significant increase of 
computing power in the last few decades makes it pos- 
sible to undertake direct numerical simulations of such 
systems. This allows to check the validity of theories in 
the most direct way. Such calculations are one of the 
subjects of this paper. 

Another, and more important point is the calculation 
based on the so-called dual variables approach (DVA), 
that provides a regular diagrammatic series constructed 
on top of the CPA. This method was originally intro- 
duced for classical lattice ensembleJ22l and correlated 
fermion^^H^ a,s an extension of the effective-medium ap- 
proximations. Recently, the DVA have been successfully 
applied for the Anderson localization problerrPS! for dis- 
ordered systems. It is worth to note that in all men- 
tioned papers an introduction of dual variables is based 
on Hubbard-Stratonovich transformation in the expres- 
sion for the partition function. For the Anderson local- 
ization problem, replica method was also utilized. In 
the present paper we show that for disordered system, at 
least at the basic level, neither the Hubbard-Stratonovich 
transformation nor replica method are not nessesary for 
the construction of the DVA. The change of variables, 
and further construction of the diagram series can be 
done directly in the expression for the linear response of 
the system. 

To highlight the effects of long-range correlations, we 
perform a simulation of a lattice gas with dipole-dipole 
interaction and compare the numarical data with results 
of the CPA, NLCPA, and the DVA. 

In the "Definitions", we describe details of the sys- 
tem under consideration and discuss methods used for 
study of the response function of those systems. In the 
next three subsections, we outline the CPA, NLCPA and 
the dual-variables approach, respectively. The "Calcula- 
tions" section is devoted to computational details, and in 
the "Results" we compare the numerical simulations and 
analytical calculations described above. 



II. DEFINITIONS 

We consider the optical response of a planar film con- 
taining spherical metallic particles. The size of particles 
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and the lattice constant are assumed negligible in com- 
parison with the wavelength of the incident radiation. 
Also for simplicity, we consider only the case of the elec- 
trical field being normal to the film plane. The general 
case can be considered similar ly^°. Particles are placed 
at the nodes of a regular square lattice (arbitrary lattice 
symmetry could be assumed), but their susceptibilities 
are random uncorrelated quantities. In the dipole ap- 
proximation, this system is described by the system of 
equations for the dipole moments df. 



d, = 



(1) 



where ai are (random) susceptibilities of the particles, 
and Cij describes their dipole-dipole interaction. Since 

we stick to the Rayleigh limit, Cij — (Kronecker 5 

ij 

is included to exclude a self-action). Note that in the op- 
tical problems susceptibility has a finite imaginary part, 
wheras in problem of the electrons localization assume an 
infinitesimal displacment of poles from the real axis. In 
this paper, we are interested in the average value of the 
dipole moment, (d). Using the ergodic hypothesis one 
can equally average over the sites of an infinite system 
or over the realizations. A distribution function of the 
susceptibility a should be specified. We assume a binary 
distribution, so that takes the value a"*" with a prob- 
ability p, and a~ with a probability (1 — p). There is an 
important special case of a~ = 0, which corresponds to 
a partially filled lattice (lattice gas) with a filling factor 
p. It also should be pointed that a has a spectral depen- 
dence a = a{uj), where w is the frequency of the incident 
field. 

The equations ([T]) can be represented in a matrix form: 



C-d) 



(2) 



In real-space, the matrix a is diagonal, with random un- 
correlated elements. The regular matrix C becomes di- 
agonal after the Fourier transform. The solution of ^ 
after averaging over system configurations gives the av- 
erage dipole moment of the whole system and can be 
represented as follows: 
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C 



(3) 



where G is the response function. The later is diagonal 
in quasimomemtum space (fc-space), due to the transla- 
tional invariance of the averaged system. Mention here 
that Re Gk=o function represents the average dipole mo- 
ment of the whole system, wheras Im Gk determines 
the density of state (DOS) of the system. Such phenom- 
ena as localization, optical scattering and high harmonic 
generation are determinded by higher order Green's func- 
tions. For a weak disorder the standard diagrammatic 
expansion is valid, in general case a more sophisticated 
consideration is needed. 



a) ^ 



= b) i 1= 




FIG. 1: Diagrammatic notations 



A. The effective-medium concept 

The most trivial effective-medium approach to calcu- 
lation of the response of a disordered system is to replace 
all the dipole moments with an average value, so that ^ 
become: 



(4) 



For the system under consideration this approach is the 
so-called virtual crystal approximation (VGA). 

The Fourier-transform gives the solution in reciprocal 
space 



qVCA 



Gk 



(5) 



One can observe that such a mean-field approximation 
corresponds to neglecting the correlations between the 
susceptibility of a particle and the acting field; math- 
ematically this is the decoupling (^cti^jCij ■ dj^ ^■ 

{ai) "^j Cij ■ (dj). Strictly speaking, this is not true even 
for a single fluctuating site placed in a regular medium. 
Indeed, the dipole moment results in a polarization of 
the medium, and that polarization contributes to the 
field acting on the site. This is similar to the effect of 
the so-called electrostatic image, that acts on a dipole 
placed near a substrate. The larger the dipole moment 
is, the larger the image field becomes. So, the suscepti- 
bility and the field become correlated, and the mean-field 
decoupling breaks down. 

Such "trivial" correlations are completely taken into 
account within the GPA. In this approximation, a hy- 
bridized single-site problem is considered, so that instead 
of (a) one calculates the average: 



9 ■ 



A 



(6) 



where 5 represents the regular self-influence of the node 
due to the environment (bath). The average response is 
now given by the Green's function 



G 



CPA 
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A-Gk 



(7) 



As a last step of the GPA formulation, one requires that 
the single-site average Q mimics the local part of the 



4 



lattice Green's function, 



(8) 



This equation allows to find self-consistently the hy- 
bridization A. 



B. The non-local CPA 

The natural way to further develop an effective 
medium theory, is the construction of a self-consistent 
scheme, which takes into account non-local correlations, 
using the self-consistency condition for a set of sites (clus- 
ter). The main disadvantage of such methods is the vi- 
olation of the translational invariance of the underlaying 
lattice as a result of the embedding of a real-space cluster 
into the medium. In the NCLPA self-consistent clusters 
with the periodical boundary conditions are constructed 
in the reciprocal space, and the effective medium in this 
method preserves the translational invariance of the ini- 
tial lattice. It is constructive to represent the NLCPA 
effective medium in coordinate space not as a set of site 
clusters, but as a set of interacting sublattices. In the 
reciprocal space this is equivalent to dividing the Bril- 
louin zone of the initial lattice into A^^^ tiles. In contrast 
with the single site approximation for the self-energy in 
the CPA, the exact self-energy in the NLCPA is approx- 
imated as a step function with the number of steps in 
reciprocal space, corresponding to the number of consid- 
ered sub-lattices: 



E(k) ~ E(K„). 



(9) 



This assumption accounts for the non-local correlations 
described by off-diagonal elements of the self-energy only 
for sites located on different sublattices. 

The next approximation made in the NLCPA is setting 
the phase factors in the Fourier transformation to unity 
for all k 7^ K„. Physically this means that each site of 
the cluster interacts with the medium irrespective of its 
position in cluster. 

The corresponding Green's function is represented by 
coarse-grained propagators (K is taken from the set 
of K„): 



5K = §EG(K + k), 



N 



(10) 



where the k - summations run over the momenta of the 
n-th cell, N is the number of nodes in the lattice, and Nc 
is the number of tiles. The explicit representation of the 
coarse-grained Green function in reciprocal space has a 
standard form: 
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N ^ (a-i) - C, 



K+k 



(11) 



It is necessary to introduce the cluster-excluded propaga- 
tor (also called the cavity Green's function) Qk to avoid 
the occurrence of a cluster self-action: 



^NLCPA 



(12) 



The self-consistent scheme in this approach is based on 
a requirement that the coarse-grained Green function 
equals the Green's function for the cluster placed in a 
corresponding cavity: 



G 



NLCPA 
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(13) 



where the average is taken over the cluster configura- 
tions. In a practical calculation, the righ-hand side of 
the expression ( 13 ) is calculated in the real space, with a 
subsequent Fouirer transform. For a small cluste and a 
binary distribution of a^, the cluster Green's function can 



be found analytically. Formulas ( TT]|l3 ) form a closed set 
of the NLCPA equations. 

To show similarity between the CPA and the NLCPA, 
one can rewrite (131 as follows: 



G 



K 



A 



(14) 



where A = gj^^ — Sk — (a 
appears to be an extension of 



In this notation ill 
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N 



E 



9k 



A-C, 



(15) 
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Clearly, the NLCPA passes into the CPA for Nc^l. In 
the limit ol Nc ^ N the results of the NLCPA become 
exact. 

The CPA contribution to the self-energy part con- 
sists of only the members corresponding to scattering 
at a single node, and consequently the self-energy is k- 
independent: = E. Because of the assumption about 
the approximation for the exact self-energy, the NLCPA 
takes into account (in an approximate way) non-local cor- 
relations at the lenghtscale comparable with a cluster 



C. Dual variables approach 

We start the discussion of the dual variables method 
by formally introducing the following expression: 



M=(C + A)-i-(a-i + A)-i, 



(16) 



where A is a diagonal, regular matrix. Inversion of this 
matrix and subsequent averaging gives = {G + 

A)(a-i - G)-\G + A) + (C + A), so that instead of 
eq. ([3| one can calculate (Af~^): 

-C)-i)(C + A) + (C + A), 

(17) 



(M-i) = (C + A)((c 
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We may assume that (M-i) can be represented as an 
expansion: 
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R + (RYR) 



(18) 



where R = {Aiy^ and Y = {g + A)-i - {{g + A)'^). 
Diagrammatic techniques can be used for representation 
and further work with this expansion. By construction, 
we require that 



{YRY) = 0, 



(19) 



at each order of the expansion (self-consistent condition). 
This condition in diagrammatic form reads (Fig. [l^), 
where the sohd hne represents the dual propagator, and 
dashed line is corresponds to influence of the disorder. 

In reciprocal space this requirement is represented as 
follows: 



E 



(Ck + A)-i-((g + A)-i) 



0. 



(20) 



At the zeroth order of the expansion (18) (Af ^) 



(M)-' and since 1" 7^ 0, it immediately follows from 
(pO| that: 



N ^ q 



(21) 



It is clear that later is the same as the self-consistent 
condition of the CPA In other words, the zeroth 

order of this approach is the single-site self-consistent 
effective medium approach, equal to the CPA. 

One can observe that the theory is indeed v ery sim- 
ilar to the dual-fermion approach^-: formulas ( |17| ) and 
( [l9| are the same as corresponding expression for the 
dual fermions, although Hubbard-Stratonovich transfor- 
mation was not used in the presented approach for a dis- 
ordered system. 

Let us now consider the leading-order correction to the 
CPA, as it appers from the expansion (18). As long as 
(Y) = 0, all odd members of this expansion vanish. The 



second term of the ( 18 ) is also equal to zero, since we re- 



quire ( 19 ) by construction of the method. It now follows 



that the next non- vanishing term of the (18) is of fourth 
order in Y. We can use standard techniques to achieve 
the following representation of this approximation in the 
form of the Dyson equation: 



G- 



R- 



(22) 



where G stands for and represents the Green 

function in dual variables; E is the dual self-energy part. 
The diagrammatic representation of this equation has the 
usual form shown on Fig. [ij;, where G is depicted as the 
double line. The self-energy part S, by the assumption 
above has the next form (Fig.[I]3). 



It is worth to mention that the CPA could be repre- 
sented as diagrammatic equal of the double line to the 
solid one. 



Taking into account the self-consistent condition (19 1, 



we can obtain an explicit form of the Green's function in 
reciprocal space: 



1 



(23) 



g-i+A-Ck-(.9 + Sk^)-i' 
where A is found using an iteration procedure. 

III. CALCULATIONS 

For the model under consideration we have performed 
calculations of the response function using different 
methods, described above. Results of the numerical sim- 
ulations are used as reference to compare with the results 
of our analytical approaches. 

Numerical simulations were performed for the square 
lattice formed by A'' = (32 x 32) nodes with periodic 
boundary conditions. The average response function is 
calculated by averaging the response functions for a set of 
random system configurations. For each system configu- 
ration the response function is calculated directly. This 
leads to the series of multiple inversion of the matrices 
with corresponding dimension N'^ , the typical number of 
realizations is 5 • 2^°. The response function for differ- 
ent system parameters can be calculated independently, 
and it is natural to perform these calculations in a par- 
allel way. Taking into account that at each realization of 
the system only the diagonal elements of the correspond- 
ing matrix are changing, a special algorithm can be used 
for the inversion of this set of matrices. This algorithm 
allows to efficiently vectorize computations and perform 
the calculations using the GPGPU (general-purpose com- 
puting on graphics processing unit). 

Calculations are organized as follows: the regular part 
of the matrix is inverted on the CPU using a standard 
numerical alghoritm, then, for each system realization 
the diagonal matrix elements are changed consequently, 
and the inversion of the matrix is performed on the GPU, 
using relation (based on the Sherman-Morrison formula) : 



G', = G, 



SM 



l + Mii- 6M 



GiiGij, 



(24) 



where G[j is inverted from matrix Mij with l-th diagonal 



element changed to 



Mn 



SM. This algorithm 



allows to perform calculations up to 70-times faster than 
corresponding computations on the single CPU. 

We calculate the dependence of the response function 
on the inverse single-particle susceptibility for systems 
with different filling fractions, and compared it with an- 
alytical results. We apply the idea of the NLCPA for 
the system under consideration with the Brillouin zone 
divided into N^, = 4: tiles ((2x2)- tiling). For systems 
with different filling fractions we performed calculations 
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of the dependence of the response function on the inverse 
susceptibiUty of a single particle, and performed compar- 
ison with analytical results. Since the susceptibility of a 
single particle is frequency dependent, these dependen- 
cies can be discussed as spectral dependencies. Let us 
discuss the results of analytical approaches applied to 
the system under consideration. 



IV. RESULTS 



For presentation purposes, it is better to switch from 
the susceptibility a to a frequency argument. We suppose 
a = ao / {uj — ujQ + 17) with a small loss factor 7. 




FIG. 2: Imaginary part of the response function for systems 
with different filling fractions. Interaction constant Co ~ 0.30. 

Thus created spectral dependence of the density of 
states (that is, local part of ImGk) is shown in Fig. [2] for 
different filling factors. One can observe that the quali- 
tative behavior of the curves is complex. The numerical 
data shows a central resonant peak and two satellites, and 
a broadening of the entire curve with an increase of the 
filling factor. As it follows from the CPA curves, the local 
single-particle physics is responsible for the broadening 
effects and the overall shape of the curve, whereas the 
satellite structure is completely missed in the CPA and 
must be therefore attributed to non-local effects. The 
displacement of the satellites from the central peak is al- 
most independent of the filling factor. It is natural to 
conclude from this observation that the satellites appear 
because of certain short-range effects, like a formation of 
resonant dimers of two particles placed at neighboring 
nodes. Indeed, long-range collective effects involving a 
large number of particles would be strongly dependent 
on the concentration. 

As it was described above, cluster methods and dual 
variables provide complimentary approaches to the prob- 
lem of non-locality. The NLCPA describes the short- 
range effects related to resonances due to some local 
configurations (for example, the formation of a resonant 
dimer of two neighboring particles, as discussed). Those 
effects are treated unperturabtively: an infinite number 
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FIG. 3: Quasimomentum dependencies of the imaginary part 
of the Green function at the different points. 



diagrams is virtually summed up. Contrary, the DVA 
allows to take into account both short- and long-range 
correlations on equal footing, but only within a certain 
perturbation procedure (in our case, just taking into ac- 
count the leading-order diagram). From this argumen- 
tation one should expect that the satellite structure is 
better described with the NLCPA, and indeed, our cal- 
culations confirm this. Similarly to the previous results 
for the Anderson model the NLCPA curves do show the 
satellites, although their positions do not always coincide 
with the numerical data. On the other hand, the DVA 
completely misses this effect. 

Now let us turn to the results for the k-dependent spec- 
tral function ImGk near the resonance of the density of 
states. Figure [3] shows the data exactly on the resonance 
and at a point on the shoulder of the curve (pointed by 
arrows in Fig. [2]). Our numerical data demonstrates a 
strong change of the dispersion law, that is not accu- 
rately captured by the CPA. Besides a quantitative dis- 
agreement, at the main resonance of the density of states 
one can observe a qualitative difference: the numerical 
data along (7r,7r) — (7r,0) line at border of the BZ show 
a much flatter spectral function than the CPA predicts. 
This flattening can be interpreted as an effect of a strong 
disorder. It occurs for certain region of the quasimomen- 
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turn values and is consequently a non-local effect. 

For the k-dependence of the spectral function, one can- 
not expect a good result from the NLCPA, because the 
(2 X 2) - tiling is too crude to describe an actual k- 
dependence, and a step-like change at the borders of tiles 
is unphysical. Indeed, one can discuss only the slope of 
the NLCPA curves in Fig. |3j the detailed behavior is not 
reproduced. It is also hard to speculate about a repro- 
duction of the flattening phenomenon. 

Contrary, the DVA describes short- and long-range 
phenomena at the same footing, and thus produces 
smooth k-dependencies. It also describes the observed 
flattening phenomenon, demonstrating a behavior quali- 
tatively very similar to the numerical data. We conclude 
that the DVA handles the physics of the strong disorder 
near the DOS peak energy. 

We should remark that such a picture is not generic 
for all points of the spectral curve. For instance, one 
should not expect a good result from the DVA at the 
satellite points, as the formation of the satellites itself is 
not handled by the DVA. Indeed, as Fig.[2]shows, neither 
the NLCPA nor the DVA are satisfactory in this respect. 
In fact, the best fit of the numerical data is given by just 
the CPA. 

More progress can be expected with use of combined 



schemes, including both cluster and dual-variables con- 
tributions. We made an effort to implement the sim- 
plest scheme of this kind. We substituted the short- 
ranged part of the self-energy calculated in the CPA by 
the corresponding part of the NLCPA self-energy. Re- 
sults of this method are shown in figures by the curves 
marked as NCLPA-I-DVA, and it is seen that this scheme 
doesn't works satisfactory: no significant improvment is 
found; the curve shows kinks and a wrong behavior at the 
shoulders. A more advanced combined method is clearly 
needed. 

We conclude that the dual- variables approach provides 
a powerfuU framework for systems with a diagonal disor- 
der. It is more suitable for the description of a long-range 
spatial correlations, whereas a short-range physics is bet- 
ter captured by the NLCPA (cluster) methods. 
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